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Thomas Ules, Daniel Liiftner, Eva Maria Reinisch, Georg Roller, Peter PuschnigQ and Michael G. Ramsejj^] 
University of Graz, Institute of Physics, NAWI Graz, Universitatsplatz 5, 8010 Graz, Austria 

(Dated: August 20, 2015) 

With angle resolved photoemission experiments and ab-initio electronic structure calculations, the 
pentacene monolayers on Ag(110) and Cu(110) are compared and contrasted allowing the molecular 
orientation and an unambiguous assignment of emissions to specific orbitals to be made. On Ag(110), 
the orbitals remain essentially isolated-molecule like, while strong substrate-enhanced dispersion and 
orbital modification are observed upon adsorption on Cu(110). We show how the photoemission 
intensity of extended systems can be simulated and that it behaves essentially like that of the isolated 
molecule modulated by the band dispersion due to intermolecular interactions. 


PACS numbers: 33.60.+q,31.15.ae,73.20.-r,68.43.-h 

I. INTRODUCTION 

The effects arising upon adsorption of 7r conjugated 
molecules on metal substrates are of great interest due 
to their importance in organic electronics. The ability to 
identify electronic states with specific molecular orbitals 
and to determine their energy ordering is vital to the un¬ 
derstanding of overlayer/substrate systems. This can be 
particularly difficult for strong chemisorptive interactions 
where significant broadening or level splitting occurs or 
when intermolecular dispersing bands are formed. In 
the last decade ab-initio electronic structure calculations, 
particulary within density functional theory (DFT), have 
become almost indispensable in the interpretation of ex¬ 
perimental results. However, it is becoming recognized 
that DFT results can be misleading owing to approxima¬ 
tions for exchange-correlation effects which may severely 
affect predicted adsorption geometries and/or the elec¬ 
tronic structure of organic/metal interfaces^. The pen- 
tacene/Cu interface is such a case. Despite numerous ex¬ 
perimental and theoretical studies, there is as yet no con¬ 
sensus on the orbital assignment. Here we demonstrate 
how the angle resolved photoemission (ARUPS) tech¬ 
nique, that is becoming known as orbital tomography®^, 
can provide a definitive assignment of the emissions even 
for strongly interacting extended two-dimensional (2D) 
systems and give insight into the nature of dispersion 
and hybridization. 

The pentacene (5A) monolayer on Cu(110) has been 
extensively studied. From their ARUPS investigations 
Seki et alF concluded that selection rule arguments can¬ 
not explain the photoemission behavior and suggested 
strong hybridization with the substrate could be mod¬ 
ifying the orbital symmetr}®. They also suggested the 
appearance of dispersing, interface-induced sta tes a ris- 
ing from substrate interaction^ 9 . Ferretti et al. i0 11 in¬ 
troduced the possibility that the mixing of 5A molecu¬ 
lar orbitals with the Cu substrate leads to ” interaction 
states localized at the interface”, where their calcula¬ 
tions suggest partial occupation of the lowest unoccupied 
molecular orbital (LUMO) whose symmetry might not 
be directly related to the original molecular state. On 
the same Cu(119) vicinal surface, scanning tunneling mi¬ 


croscopy (STM) images also suggest partial occupation 
of the LUMO and dispersive electronic states associated 
with a perturbed electron charge density distribution. 12 
Also a recent combined DFT and ARUPS study of the 
5A/Cu(110) system concludes that only partial LUMO 
occupation takes place! 13 . 

The simple relationship between the angular distribu¬ 
tion of the photoemission current and the Fourier trans¬ 
form of the emitting molecular orbital has been shown to 
be reasonable for a number of mol ecular ad sorbate sys¬ 
tems on various noble metal surface ffl^- 9 ^^! This allows 
molecular orientations to be determinecP^, molecular 
orbital energy ordering to be deduced 1 ^ and even the 
reconstruction of the molecular orbitals in real spacd^. 
For these systems, the angular-dependent emissions can 
essentially be accounted for by the photoemission from 
isolated molecules, thus intermolecular orbital overlap 
plays a minor role. However, in extended 2D overlayer 
systems, a description in terms of isolated molecular or¬ 
bitals is no longer strictly appropriate. Here, with the 
comparison between pentacene monolayers on Ag(110) 
and Cu(110), we show how orbital tomography and the 
Fourier transform description can be applied to extended 
systems with strongly dispersing emissions. In so doing, 
we provide a definitive description of the 5A/Cu ARUPS 
and show that the LUMO is in fact fully occupied and 
displays a substrate-induced dispersion which is signifi¬ 
cantly larger than that reported for similar organic over¬ 
layer systemd^ 17 . Moreover, unlike all previous tomog¬ 
raphy studies we here show evidence for a modification 
of orbital shape on adsorption. 


II. EXPERIMENTAL DETAILS 

Photoemission experiments were performed at 
BESSYII using a toroidal electron-energy analyzer 
described previousl}^ which was attached to the 
beamline U125/2-SGM of the synchrotron radiation 
facility BESSY II, Helmholtz-Zentrum-Berlin. Photon 
energies of 30 and 35 eV and an incidence angle of x = 
40° with respect to the surface normal were used. The 
polarization direction is in the specular plane, which 
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also contains the photoelectron trajectory measured. 
Emitted photoelectrons are recorded simultaneously 
with polar angles 0 of —80° to +80° with respect 
to the surface normal in an energy window of 1 eV. 
The energy window is divided into 40 individual slices 
which goes well below the analyzers intrinsic energy 
resolution of 150 meV. Note, however, for the presented 
momentum maps we take only the data from the side of 
the emission direction that is pointing in the orientation 
of the electric field vector of the incident photons, this 
is referred to as the positive side (see Fig. [l]), i.e., 6 
= 0° to +80°. This maximizes the polarization factor 
appearing in the photocurrent cross section (see Eq. [4| . 
It is noted that the molecular features are enhanced 
relative to the substrate emissions on the positive side 
(compare Fig. Hp). To obtain the full ( k^, k y ) range 
for the presented momentum maps at constant binding 
energy, azimuthal scans are made by rotating the sample 
around the surface normal in 1° steps for an azimuthal 
angle range > 180° and then imposing the substrates 
twofold symmetry to obtain the full 360°. The angular 
emission data are then converted to parallel momentum 
components k^ and k y using the relation 
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to create the momentum maps. The Cu(110) and 
Ag(110) substrates were prepared in the conventional 
way by a sequence of sputter-annealing cycles. 5A 
molecules were evaporated from an effusion cell onto the 
surfaces at room temperature with the amount moni¬ 
tored by a microbalance. The resulting monolayer LEED 
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for pentacene on Cu(110). In both cases 
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mild annealing improved the order. On Cu(110) the 
monolayer was annealed at 200° C which is above the 
5A sublimation temperature whereas on Ag(110) the 
molecule substrate bond is weaker and the annealing tem¬ 
perature must not exceed 140° C so as not to desorb from 
the monolayer. 


III. COMPUTATIONAL DETAILS 
A. Density functional calculations 

All theoretical results presented here are obtained 
within the framework of density functional theory (DFT) 
using the VASP code 19 20 . We have performed three 
types of calculations: firstly for the isolated pentacene 



FIG. 1. Sketch of an ARUPS experiment declaring the ge¬ 
ometry of the experimental setup. The incoming photon with 
the energy hv, the incidence angle x an d vector potential A 
excites an electron from the initial state ipi to the final state 
ipf. This final state is characterized by the kinetic energy Eki n 
and the momentum vector k and the outgoing photoelectron 
is detected as a function of Eki n and emission direction, de¬ 
fined by the polar angle 0 and the azimuthal angle (/). Where 
forward emissions are denoted for by a positive parallel mo¬ 
mentum and backward by a negative, respectively. 


molecule, secondly for a two-dimensional, extended free¬ 
standing layer of pentacene, and thirdly for monolayers 
of pentacene adsorbed on Ag(110) and Cu(110) surfaces. 

The isolated molecule calculations were performed us¬ 
ing a supercell with a minimum of 15 A vacuum be¬ 
tween pentacene’s periodic replica. We use the gener¬ 
alized gradient approximation (GGA) 21 I for exchange- 
correlation effects, and the projector augmented wave 
(PAW) method 22 . The simulated momentum maps of the 
pentacene HOMO and LUMO orbitals shown in Fig.[2]are 
obtained as Fourier transforms of the respective Kohn- 
Sham orbitals as described previously—. 

Electronic structure calculations for the freestanding 
monolayer of 5A have been carried out using the repeated 
slab approach with a vacuum layer of 20 A between ad¬ 
jacent layers. The generalized gradient approximation 
(GGA) 21 is used for exchange-correlation effects, and 
the projector augmented waves (PAW)P approach was 
used allowing for a relatively low kinetic energy cut-off 
of about 400 eV. We use a Monkhorst-Pack 6x12x1 grid 
of /c-points 23 ,_and a first-order Methfessel-Paxton smear¬ 
ing of 0.1 eV—. The simulation of momentum maps for 
extended systems will be described in the subsequent sec¬ 
tion. 

Finally, we have also performed DFT calculations of 
pentacene monolayers adsorbed on Ag(110) and Cu(110). 
The substrate is taken into account within the repeated 
slab approach by using five metallic layers with an addi¬ 
tional vacuum layer of 15 A between slabs. To avoid spu¬ 
rious electrical fields, a dipole layer is inserted in the vac¬ 
uum region 25 . In the case of 5A/Ag(110), we have taken 
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into account the experimental LEED structure 

mentioned above and relaxed the atomic positions of the 
molecule and the first metallic layer considering van-der- 
Waals interactions by employing the empirical correction 
scheme according to Grimme 26 . Exchan ge c orrelation ef¬ 
fects were treated either within the GGA 2 0or within the 
hybrid functional HSE 27 lwith k-point meshes of 9 x 6 x 1 
and 6 x 4 x 1, respectively, and a first-order Methfessel- 
Paxton smearing of 0.1 eV. For 5A/Cu(110), we have 



chosen a commensurate structure 


6 0 
0 2 


and an ad¬ 


sorption site similar to a previous stud}/ 13 . Also, here 
the electronic structure is calculated within the GGA and 
HSE using k-point meshes of 9 x 12 x 1 and 4x6x1, 
respectively, and using a first-order Methfessel-Paxton 
smearing of 0.1 eV. 


ulus of the Fourier transform of the initial state wave 
function ^ nq modulated by the weakly angle-dependent 
factor | A • k| 2 which depends on the angle between the 
polarization vector A of the incoming photon and the 
direction k of the emitted electron. 

Note Eq. [4] can be applied to single molecules as well 
as to extended states such as organic layer adsorbed on 
metallic surfaces. In the former case the summation of 
the Brillouin zone (BZ) reduces to just one point, the T 
point, while for the latter situation, the possible dispers¬ 
ing bands are taken into account by the band structure 
£ nq and Bloch states ^ nq and an appropriate sampling 
of the Brillouin zone. 


IV. RESULTS 
A. Momentum maps 


B. Simulation of ARUPS maps 

The Kohn-Sham energies £ nq and orbitals ^ nq of the 
relaxed structures serve as input for the subsequent sim¬ 
ulation of ARPES intensity maps within the one-step 
model of photoemission 28 . Here, the angle-resolved pho¬ 
toemission intensity 1(6, (j); Eki n , cj) is a function of the 
azimuthal and polar angles 0 and <fi, respectively, the ki¬ 
netic energy of the emitted electron E^[ n , and the energy 
ljo and polarization A of the incoming photon: 

occ BZ 

1(9,<j>, £ kin ; w, A) « ]T ]T \(ip f (6, 0; E kin ) |A • pIVWI" 

n q 

X ^(^nq ± ± (3) 

This formula can be viewed as a Fermi’s golden rule ex¬ 
pression, in which the photocurrent I results from a sum¬ 
mation over all occupied initial states ^ nq , characterized 
by the band index n and Bloch vector q, weighted by 
the transition probability between the initial state and a 
final state. For the transition operator A • p, the dipole 
approximation is assumed, where p and A are the mo¬ 
mentum operator and the vector potential connected to 
the incoming photon. The 5 function ensures energy con¬ 
servation where <f> denotes the work function. 

We further approximate the final state ijjf by a plane 
wave 29 . As outlined in more detail in a previous papei^, 
and also noted earhei^MikU, this approximation allows us 
to greatly simplify the evaluation of the matrix element 
appearing in Eq. [3j If we denote the wave vector of the 
of the final, free-electron state by k, thus E^ n = A_/c 2 , 
Eq. [ 3 ] simplifies to 

occ BZ 

ky> Ekirn ^i A) ~ EEiA-ki 2 |<e-ivwf 

n q 

x $(^nq ± 4> ± Ekin — &). (4) 

We obtain the simple result that the matrix element from 
a given initial state nq is proportional to the square mod¬ 


Monolayers of 5A on Ag(110) and Cu(110) were char¬ 
acterized by low energy electron diffraction (LEED) and 
ARUPS. Fig. [2] shows LEED images of the 5A monolay¬ 
ers on the two substrates Ag(110) and Cu(110) as well 
as the momentum maps taken at the binding energies of 
the prominent adsorption induced emissions at 0.15 and 
1.2 eV on Ag and 0.8 and 1.5 eV on Cu, respectively. 
These maps clearly display the character of the LUMO 
and HOMO of pentacene as can be seen from the com¬ 
parison between the measured (exp) and the correspond¬ 
ing simulated maps computed for the isolated molecule 
(sim). The simulated LUMO map, for instance, is char¬ 
acterized by a major lobe at (0,1.4) A -1 and minor lobes 
at (±1.8, 1.2) A -1 . Indeed the data on Ag has recently 
been used to reconstruct these orbitals in real space in ex¬ 
cellent agreement with calculated orbitals for the isolated 
molecule 6 . The maps also unambiguously reveal the ori¬ 
entation of the molecules: on Ag(110) the data implies 
flat lying 5A oriented parallel to the [001] azimuth, while 
on Cu the molecules orient along the [1-10] azimuth as 
depicted in Fig. and h, respectively. It should be 
noted that the structure and density of the two mono- 
layers are very similar, however, on Ag the molecules lie 
across while on Cu they he parallel to the close packed 
rows of the substrate. 

Thus, the ARUPS momentum maps of Fig. [^immedi¬ 
ately clarify the 5A/Cu system. In contrast to specula¬ 
tions in the literature, the symmetry of the orbitals is not 
modified significantly on hybridization. This allows the 
emissions to be unambiguously assigned. The features at 
0.8 and 1.5 eV are the fully occupied LUMO and HOMO 
and not the HOMO and HOMO-1 emissions as previously 
and, in the light of our results, erroneously assigned on 
the basis of a comparison between DFT-calculated and 
experimental binding energies 13 . On closer inspection 
of the LUMO map on Cu, however, an internal struc¬ 
ture and a change in the ^-position of the minor lobes 
becomes visible that will be argued to be the result of 
intermolecular dispersion and changes in the orbital size, 
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FIG. 2. Panels a) and e) show LEED images of 5A mono- 
layers on Ag(110) and Cu(110), respectively, b) and c) show 
ARUPS momentum maps of 5A/Ag(110) at binding energies 
0.15 and 1.20 eV (exp), respectively, compared to simulated 
LUMO and HOMO maps of the isolated 5A molecule (sim). 
Panels f) and g) show the corresponding data for 5A/Cu(110) 
at binding energies of 0.80 and 1.50 eV. Panels d) and h) 
show structural models of 5A/Ag(110) and 5A/Cu(110), re¬ 
spectively, as deduced from LEED and the momentum maps. 
In h) the orbital structure of the LUMO is overlaid. 


respectively. 


B. DFT results for adsorbed monolayers 

Before analyzing this internal structure of the LUMO 
in more detail, we discuss DFT results of the pentacene 
monolayers on the Ag(110) and Cu(110) surfaces in com¬ 
parison with ARPES data. Starting from the relaxed 
structure of 5A/Ag(110), we have computed the density 
of states projected onto the orbitals of a freestanding 
pentacene layer, both, within GGA and self-consistently 


pentacene / Ag(110) 



E - E f (eV) 


pentacene / Cu(110) 



E - E f (eV) 


FIG. 3. Density of states of a monolayer of pentacene on 
Ag(110) projected onto orbitals of the free pentacene molecule 
using PBE-GGA (a) and the hybrid functional HSE06 (b). 
Projections onto the LUMO (blue), the HOMO (red), the 
HOMO-1 (orange) and the HOMO-2 (green) are shown. For 
comparison, panel (c) shows experimental ARPES data along 
the [001] azimuth (black solid line) and an azimuthal direc¬ 
tion 45° between [001] and [1-10]. The colored curves are 
obtained by /^-integrating experimental momentum maps 
identifying the orbital character of the emissions. Panels (d)- 
(f) show the corresponding data for a monolayer of pentacene 
on Cu(110). 


within the hybrid functional HSE 27 32 . It is known that 
the incorporation of a fraction of exact exchange in the 
hybrid HSE mediates self-interaction errors thereby im¬ 
proving the orbital energies and thus the description of 
adsorbate system 35 . The results are shown in pan¬ 
els (a) and (b) of Fig. [3j For comparison, panel (c) 
shows experimental ARPES data recorded over a large 
energy window along the two different azimuths [001] 
(black solid line) and in a direction 45° between [001] 
and [1-10] (black dashed line). Panel (c) also includes 
the k x /^-integrated ARPES data of the momentum maps 
that have been recorded over four energy windows of 
1 eV centered around the molecular emissions of the 
LUMO, HOMO, HOMO-1, and HOMO-2, respectively. 
When comparing GGA with HSE results, we observe 
that the HOMO, HOMO-1 and HOMO-2 features are 
shifted to somewhat larger binding energies when using 
the hybrid functional while the position of the LUMO or¬ 
bital slightly below the Fermi level remains unchanged. 
Compared to experiment, we note that the HSE calcu¬ 
lation clearly improves energy position of the HOMO-1 
and HOMO-2 and also leads to a slightly better agree¬ 
ment of the HOMO position. Regarding the LUMO, both 
functionals GGA and HSE indicate partial occupation in 
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good agreement with experiment. 

For the case of Cu(110) (Fig. [3j panel d-f), HSE again 
yields much better agreement with experiment. However, 
the LUMO, of prime importance in this study, is still in 
poor agreement with experiment. Its computed position, 
both, obtained within GGA and within HSE, is in error 
compared to experiment as the LUMO is located at the 
Fermi level and not fully occupied with only a slight im¬ 
provement gained by the hybrid functional calculation!^. 
Note that various commensurate structures and adsorp¬ 
tion sites did not change this finding. On the other hand, 
HSE does considerably improve the energy position of 
the deeper lying orbitals, for instance, the HOMO. As 
a side-note: It was exactly the wrong energy position 
of the HOMO in the GGA result which had motivated 
Muller et al. to erroneously assign the LUMO emission 
to the HOMO. Note that the computed DOS of their 
paper 1 ^ agrees with our GGA-DOS and their experimen¬ 
tal ARUPS data is in line with our experiment. 


C. ID Band dispersion 

Let us now concentrate on the energy dispersion of 
the LUMO state as suggested by the internal structure 
of the momentum maps. Fig. [4] shows the E vs. k re¬ 
lation (bandmap) along the long molecular axes of 5A 
for a) a monolayer an Ag(110), b) a half monolayer on 
Cu(110) and c) a complete monolayer on Cu(110) from 
the Fermi edge (Ep) down to 2.2 eV binding energy. On 
Ag, the bandmap shows a feature located at the Fermi 
edge at a k value around 1.35 A -1 whose intensity grad¬ 
ually tales off with increasing binding energy indicating 
that the LUMO is half-filled. Note that in this energy 
range, there are no apparent changes in the momentum 
maps. In Fig. EL the half monolayer of 5A on Cu(110) 
is investigated. By examining such a low coverage, with 
no long-range order apparent from LEED, any effects of 
intermolecular dispersion are minimized and only effects 
due to the molecule-substrate interaction are expected. 
The bandmap of the 5A submonolayer on Cu(110) in 
the [1-10] direction (Fig. 2b), i.e. along the long molec¬ 
ular axes, shows, besides the Cu d and sp bands, emis¬ 
sions from the major 5A LUMO lobe visible around 1 eV 
binding energy. Comparing with Ag, the LUMO shifts 
from Ep and being half-occupied to a binding energy of 
roughly 1 eV and being fully occupied on Cu. This is 
indicative of the strong bonding interaction of 5A with 
Cu, nevertheless, the band map shows no obvious sign of 
intermolecular dispersion. 

Upon the formation of the well ordered monolayer on 
Cu, the situation is changed with new structure appear¬ 
ing in both the band and momentum maps. As evident 
in Fig. H, the principle molecular induced emission now 
oscillates in energy from Ep down to 1.2 eV in the range 
of fc[i_io] from 1.0 to 1.9 A -1 suggestive of strong inter¬ 
molecular dispersion (see Fig. [5 Jd for an enlarged view). 
To emphasize the E vs. k relation of the LUMO along 



FIG. 4. a) Bandmap of a 5A monolayer on Ag(110) in the 
[001] azimuth, b) and c) show band maps in the [1-10] az¬ 
imuth of a half and a complete monolayer of 5A on Cu(110), 
respectively. Vertical lines indicate the Brillouin zone bound¬ 
aries (yellow) and the F points (purple) of the 5A overlayer. 

the long molecular axes direction, the Brillouin zone (BZ) 
boundaries and F points of the 5A monolayer structure 
deduced from LEED are indicated in Fig. H- Note that 
without differentiation or other enhancements of the raw 
data, the large LUMO dispersion with the periodicity of 
the overlayer is clearly evident in the 4th and 5th zones. 
The high binding energy side of the band is seen to be 
at T and the low binding energy at the zone boundaries. 
Given the LUMO orbital’s topology (Fig. [2ja) , the band 
dispersion running up from F to the BZ boundary is qual¬ 
itatively consistent with a chain of pentacene molecules. 

D. 2D Band dispersion 

Having qualitatively understood the one-dimensional 
dispersion along the long molecular axes, a more quan¬ 
titative and, above all, an understanding of the two- 
dimensional (2D) LUMO dispersion on Cu and its mani¬ 
festation in the photoemission experiment is desired. To 
this end, we first discuss the computed 2D LUMO disper¬ 
sion of the free-standing layer which is plotted in Fig. 
with the band energy E(k^ 0 i] > &[i-io]) color-coded: green 
indicating the center of the band, and red (blue) the top 
(bottom) of the band, respectively, where for reasons of 
clarity just one of the mirror domains is included. Note 
that although our 5A/Cu(110) HSE calculation includ¬ 
ing the substrate considerably improves the agreement 
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with experiment over GGA results, we refrained from 
using it here for the analysis of the experimentally ob¬ 
served band dispersion of the LUMO for the following 
reason. The experimental overlayer structure is non- 
commensurate while, in order to enable a DFT calcu¬ 
lation including Cu(110), a commensurate surface unit 
cell had to be imposed (similar to a previous stud}) 13 -!). 
This difference would result in distinctly different period¬ 
icities in the measured and simulated momentum maps 
preventing a one-to-one comparison. In contrast, the 
free-standing layer simulation with the correct structure 
enables a one-to-one analysis of the experimentally ob¬ 
served LUMO dispersion. 



FIG. 5. a) Color-coded two-dimensional dispersion 

F(/c[ooi], fc[i-io]) °f the LUMO for a freestanding 5A mono- 
layer with the structure of Fig. g}i. Superimposed in white 
are the Brillouin zones of Cu(110) (rectangle) and the 5A 
monolayer (rhombus), and in black an intensity iso-line of the 
calculated isolated 5A LUMO momentum distribution, b) 
Close-up of the dispersing LUMO band with red, green and 
blue lines indicating the energy positions at which momentum 
maps are extracted for comparison with the calculated pho¬ 
toemission momentum distributions in c) top d) middle and 
e) bottom of the band. Note that in panels c)-e) simulated 
momentum maps ( sim) are compared to experimental maps 
[exp) as detailed in the text. 


Inspection of Fig. Hk shows that the character of the 
dispersion, i.e., the /^-positions of the tops and bottoms 
of the band, may well be understood by knowing (i) the 
respective Brillouin zone, and (ii) the nodal structure of 
the respective molecular orbital in conjunction with the 
inter molecular arrangement of molecules. The first de¬ 
termines the reciprocal periodicities and the second the 
direction of the dispersion, i.e., whether the band would 
run up or down from T to the BZ boundary. The main 
dispersion is along the long real space unit cell vector, 
z.e., roughly along the long molecular axes. The calcu¬ 
lated energy spread of 20 meV for the free-standing layer 
is, however, very small and more than an order of mag¬ 
nitude lower than the experimentally observed one. 

In the ARUPS experiment, one does not observe the 
2D dispersion of the LUMO-derived band over the en¬ 
tire /c-range due to matrix element effects. For isolated 
molecules, the selection rules are well-described by the 
FT of the isolated molecular orbitals. In analogy, in or¬ 
der to simulate the photoemission intensity distribution 
in a quantitative manner, we thus need to consider the 
wave functions of the extended 2D molecular system and 
compute their Fourier transforms according to the plane- 
wave final state approximation 1 ^. Evaluating Eq. [ 4 ] for a 
free-standing layer, we simulate momentum maps (sim) 
at the top, middle and bottom of LUMO-derived band 
and compare them with the experimental maps (exp) 
taken at the 0.15 eV, 0.60 eV and 1.0 eV binding en¬ 
ergy in panels c, d and e of Fig. |5j respectively. These 
energy positions are also indicated as horizontal lines in 
the enlarged band map shown in panel b. Note that 
in the simulations, we have taken into account that the 
monolayer structure consists of two mirror domains. 

The experimental and simulated momentum maps are 
in remarkably good agreement. For the major lobe, in 
going from top to bottom of the band, not only is the gen¬ 
eral shape in agreement, one also observes an increase in 
the extension of the feature in fc[ 0 01] and a splitting of the 
emission in /c[!_ 10 ]. Although weak there is also agree¬ 
ment in the behavior of the minor lobes. For instance, 
the shift in /c[!_ 10 ] from 1.0 to 1.3 A -1 is seen in both 
experiment and simulation. Also, in the middle of the 
band the splitting of the minor lobe seen in the experi¬ 
ment, although difficult to discern, is also observable in 
the calculated map. We conclude that the photoemis¬ 
sion intensity of the extended system is essentially that 
of the isolated molecule modulated by the intermolecular 
dispersion. This naturally implies that the practice of 
searching for dispersion in extended systems outside the 
k -range expected for the isolated molecules is question¬ 
able. This is illustrated in the calculated 2D dispersion of 
Fig# by overlaying the isoline of the computed ARUPS 
map of the the isolated 5A LUMO. 

While giving extremely good agreement in the three 
momentum maps shown in Fig. Ii -e, the exact E vs. 
k relation is naturally problematic when analyzed in 
terms of a free-standing layer. Indeed, as experience 
with other systems has shown, the substrate may or 
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may not (depending of the relative energy position of 
molecular and substrate levels) enhance the dispersion 
as e.g. for PTCDA/Ag( 110 )P^ and NTCDA/Ag(110) or 
NTCDA/Cu( 100 ) 17 . But in all these cases the charac¬ 
ter of the dispersion is left unaltered and therefore the 
appearance of momentum maps is also unchanged when 
going from the free-standing layer to the adsorbed mono- 
layer. The magnitude of the dispersion is generally in¬ 
fluenced by the molecule-substrate interaction. The case 
of pentacene/Cu( 110 ) stands out because the observed 
dispersion is the largest one measured so far and may 
thus be even termed substrate- induced since the band 
width of the free-standing layer’s LUMO band is almost 
negligible. 

Indeed, our calculation for 5A/Cu(110) using the com¬ 
mensurate structure does show enhanced dispersion com¬ 
pared to the free-standing layer without the substrate, 
however, not as large as the experimental observation 
presumably as a result of the underestimation of the 
LUMO binding energy. 


E. Molecule-Metal Hybridization 

For the 5A/Cu system, in addition to the structure in¬ 
troduced by dispersion, there is also a distinct change in 
the ^-position of the minor lobes in the /q 0 oi] direction. 
In Fig. [ 6 j simulated momentum maps of the LUMO for 
the gas-phase and for the 5A/Cu(110) adsorbate system 
are compared to experimental maps of the 5A LUMO on 
both surfaces. On inspection of Fig. [ 6 ]i , one sees that 
on Cu, the minor lobes (green crosses) have maxima at 
k y = ±1.5 A -1 which is significantly lower than, both, 
the simulated values for the isolated molecule (b) and 
the experimental value on the Ag surface (e) (± 1.8 A -1 ). 
This shift is a consequence of the molecule-copper inter¬ 
action and can be interpreted as a ~ 20 % increase of 
the lateral orbital size upon adsorption. It can be ratio¬ 
nalized by recalling the reciprocal relationships between 
the LUMO shape in real space and the corresponding 
momentum space patterns as illustrated in Figs. [ 6 ^i and 
b. Here the maxima of the major and minor lobes are 
marked as red and green crosses in the momentum map, 
respectively, while the corresponding real space dimen¬ 
sions are indicated by scale bars of length 2ir/k. We 
note that the width of the LUMO orbital along the short 
molecular axis is reflected by the k y maximum of the mi¬ 
nor lobe. Thus, unlike all molecular adsorbate systems 
so far reported 4 * 6 * 14 * 37 *, the orbital tomography of the 5A 
LUMO on Cu( 110 ) shows a distinct modification of the 
orbital’s shape from the isolated molecule. 

Our interpretation of the experimental data in terms 
of a spatial distortion of the orbital is further supported 
by our DFT calculations. Starting from the orbital en¬ 
ergies and wave functions of the HSE calculation for 
5A/Cu(110), we simulate an ARPES intensity map of the 
LUMO as shown in Fig. [ 6 J 3 . These simulations indeed in¬ 
dicate a shift of the minor lobe peak (green cross) to lower 



FIG. 6. a) Calculated LUMO of an isolated 5A molecule 
in real space, scale bars mark characteristic dimensions as 
detailed in the text, b) and c) are simulated momentum maps 
for a free 5A molecule and the monolayer 5A/Cu(110), while 
d) and e) show experimental ARPES maps of 5A/Cu(110) and 
5A/Ag(110), respectively. The red and green crosses indicate 
the position of the maxima of the LUMO’s major and minor 
lobe, the horizontal green line is a guide for the eye. 


k -values upon adsorption on Cu( 110 ). We note, however, 
that this change is smaller than the experimental obser¬ 
vation. The discrepancy is likely due to a shortcoming 
of the DFT results regarding the LUMO binding energy 
(too low) and width (too narrow) at the HSE-DFT level 
indicating an underestimation of the interaction strength 
of the LUMO with the substrate. The pronounced differ¬ 
ence between the minor lobe peak position of 5A/Cu(110) 
compared to 5A/Ag(110) can thus be seen as direct ev¬ 
idence for a strong molecule-metal hybridization for the 
former. 

While hybridization is generally observed by indirect 
means in effects on molecular emissions such as energy 
broadening of orbital emissions 16 or reduced emission in¬ 
tensity in gaps of bulk stated, we here provide direct 
evidence for changes in orbital shape. Moreover, further 
signatures of hybridization can be observed as additional 
Cu-sp-like features appearing with the overlayer period¬ 
icity (Fig. §• These are best seen on the negative k 
side where the molecular features are weak as shown in 
Fig#. These emissions are not simply the bulk sp bands 
scattered by the overlayer since mere scattering would 
replicate the Cu sp band over the entire energy range. 
Instead one observes that they do not extend up to Ep 
but stop at 1.2 eV binding energy, just below the onset of 
the LUMO orbital. We therefore suggest they are due to 
interfacial Cu sp bands hybridizing with the 5A LUMO 
such that states appearing in these bands in the 0-1 eV 
binding energy range have changed their character from 
























Cu sp to that of the molecular LUMO. 

V. CONCLUSION 

In summary, a very strong substrate-induced disper¬ 
sion of the LUMO orbital of 5A/Cu(110) along the Cu- 
row direction is traced back to a significant hybridiza¬ 
tion between the organic and the metallic states, while 
the LUMO of 5A/Ag(110) exhibits only a minor inter- 
molecular dispersion. By making use of the reciprocal 
relation between the structure of real space orbitals and 
features in the momentum maps, we are able to deduce a 
significant geometrical modification of the LUMO orbital 
upon adsorption on Cu(110), while it remains essentially 
free-molecule like on Ag(110). We believe that the main 
difference between the adsorption behavior of pentacene 
on these two surfaces arises from the distinct molecular 
orientation with respect to the close-packed metal rows. 
While on Ag(110) pentacene orients perpendicular to the 
rows, it aligns parallel to the metal rows on Cu(110). Be¬ 
cause the electronic structure, e.g. the band dispersion of 
the metallic sp-bands or the location of the surface state 
in the surface Brillouin zone, is distinctly different along 


the rows [110] and perpendicular to them [001], also the 
hybridization with adsorbed elongated species can be ex¬ 
pected to be different. 

In conclusion, for the case of pentacene monolayers 
on Ag(110) and Cu(110) surfaces, we have demonstrated 
the power of the orbital-tomography method using ex¬ 
tensive angle-resolved photoemission data for revealing 
the electronic structure of such two-dimensional organic- 
metal interfacial layers. By generalizing the theoretical 
description of the photoemission process from isolated 
molecule systems to extended two-dimensional systems, 
we are able to explain the ” fine-structure” in the exper¬ 
imental momentum maps and obtain a comprehensive 
understanding of the band dispersion of organic-metal 
inter facial layers. 
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